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Abstract 

We develop for dipole-forbidden transition a dynamical theory of two-photon paired superra- 
diance, or PSR for short. This is a cooperative process characterized by two photons back to 
back emitted with equal energies. By irradiation of trigger laser from two target ends, with its 
frequency tuned at the half energy between two levels, a macroscopically coherent state of medium 
and fields dynamically emerges as time evolves and large signal of amplified output occurs with 
a time delay. The basic semi-classical equations in 1+1 spacetime dimensions are derived for the 
field plus medium system to describe the spacetime evolution of the entire system, and numerically 
solved to demonstrate existence of both explosive and weak PSR phenomena in the presence of 
relaxation terms. The explosive PSR event terminates accompanying a sudden release of most en- 
ergy stored in the target. Our numerical simulations are performed using a vibrational transition 
X^Ti'^v = 1 —7- of para-H2 molecule, and taking many different excited atom number densities 
and different initial coherences between the metastable and the ground states. In an example of 
number density close to O[10^^]cm~^ and of high initial coherence, the explosive event terminates 
at several nano seconds after the trigger irradiation, when the phase relaxation time of > O[10] 
ns is taken. After PSR events the system is expected to follow a steady state solution which is 
obtained by analytic means, and is made of many objects of field condensates endowed with a 
topological stability. 
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I. INTRODUCTION 

Since the early suggestion [l| a variety of coherent two-photon processes have attracted 
much interest, both from theoretical |2|, J3[, [4] and experimental sides |5|, 6|, [7|, 8|, 9|. 
Our present work is focused on a different aspect of coherent two photon emission from 
A— type three level atoms (or molecules) where transition between two lower levels is dipole- 
forbidden (see Fig([T]) for the level structure). As pointed out in |10J], a macroscopic target 
made of metastable atoms in the |e) of Fig([T]) may induce a characteristic event of macro- 
coherent two photon emission, two photons exactly back to back emitted with equal energies. 
We use for this phenomenon the terminology of two-photon paired superradiance, or PSR 
in short. The term paired is used because two emitted photons are highly correlated in 
their momenta and spin orientations (most clearly seen in J = — )■ transition). The rate 
enhancement factor in the momentum configuration of the back to back emission is expected 



much larger than in the usual superradiance (SR) case llj due to lack of the wavelength 
limitation there: the coherent volume for SR is limited with the wavelength A by \^L where 
L is the target length for a cylindrical configuration, while the macro- coherent PSR has the 
coherent volume of entire cylinder irradiated by trigger. 
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FIG. 1: (Color online) A— type atomic level for PSR. Dipole forbidden transition |e) — )• I5) +7 + 7 
may occur via strong El couplings to \p). 

The usual single photon superradiance occurs irrespective of absence or presence of trigger 
due to the intrinsic instability of exponential spontaneous decay caused by dipole- allowed 
transition. On the other hand, two-photon emission occurs with a much smaller rate in 
higher order of perturbation beyond dipole-forbidden transition, hence the use of trigger 



is essential to assist the macro-coherence development for two-photon process and induce 



rapid PSR events of large signal. Quantum initiation such as proposed for SR in 12|] is 
not needed, since PSR is more akin to the triggered SR llSi], which makes appropriate the 
following semi-classical treatment. ^ 

A bas,c formalism of two-photo,> process already ex,sts. [3] for propagafou equation arrd 
[14i | for analytic results of the propagation problem and PSR emission treated as pertur- 
bation. But this formulation turns out insufficient to dynamically discuss two (back to 
back) mode propagation incorporating PSR, which seems essential for dense medium. In 
the present work we shall be able to derive a fundamental set of semi-classical equations for 
the two mode and further present formulation of two color problem as well. The essential 
ingredient in our work is derivation of a more general quantum mechanical equation both for 
the medium (Bloch equation) and the electromagnetic field (extension of quantum Maxwell 
equation to include two photon process). Only after elucidating the nature of quantum state 
of fields and medium, namely, time evolving electric field condensate, we shall go on to the 
semi-classical equation. This way we determine how two back to back modes are precisely 
coupled beyond the perturbation theory. 

We ignore the granularity and rely on the continuum limit formulation of atom distribu- 
tion, taking one spatial dimension alone, because the whole event is highly focused on one 
direction of irradiated trigger field taken as x direction. The system of semi-classical partial 
differential equations thus derived is highly non-linear, and must generally be analyzed by 
numerical simulations. This way we find explosive and weak PSR phenomena and under 
what conditions these may occur. 

Despite of its complicated non-linearity the system allows soliton solution of two kinds, 
which is obtained as steady state solutions of this non-linear system of fields and medium. 
Solitons here, in their field part, are electric field condensate which may or may not be 
moving: there can be static field condensate. The stability (against two-photon emission) 
of solitons is ensured by a topological quantum number as explained in the text below. Our 
conjecture, which is supported by numerical simulations, but not estabhshed by a more 
rigorous method, is that field condensate formed after rapid PSR phenomena is made of 
many topological solitons. After formation of field condensates, namely a stable target state 
against two photon emission, the light may propagate almost freely. The condensate state of 
field plus medium thus formed may be very useful to detect a much weaker process such as 



radiative neutrino pair emission (RNPE) [l5|, because the condensate is not stable against 
RNPE. 

A related propagation and soliton formation problem in the single photon case is the 



phenomenon of self-induced transparency (SIT) [16|] and electromagnetically induced trans- 
parency (EIT) [8i] presumably related to solitons of the kind of 17|, both at a resonant 
frequency. Both of these transparency phenomena thus appear to be directly related to 
formation of stable solitons of different kinds from ours. 

For numerical computations below, we use parameters relevant to a good target candidate 
for PSR detection, para-H2 molecule. We have in mind using para-H2 vibrational transition 
of X^S^f = 1 — 7- (X being the electronically ground molecular state). Many other atoms 
and molecules are conceivable for PSR experiments. The characteristic length scale for 
large effects is ~ 14cm, and the time scale ~ 0.5 ns for para-H2 of a molecule density of 
n = 10^°cm~'^. The number density dependence of these characteristic parameters is oc 1/n. 
We include relaxation effects of two time constants in the range of T2 > 10 ns (a feasible 
value experimentally) and Ti ^ T2 in our analysis. Origin of these relaxation constants is 
left unexplained, and this way one may use values experimentally measured by other means. 
We perform extensive numerical simulation in order to clarify experimentally observable 
PSR signals and condensate formation in forthcoming experiments. It is demonstrated that 
explosive PSR emission occurs for long targets even by weak trigger when initial coherence 
between states, |e) and \g), is present. We have identified two different types of PSR events 
caused by trigger irradiation: (1) explosive PSR in which most of the stored energy in the 
initial metastable state |e) is released as a short pulse of some time structure, and (2) weak 
PSR in which the output energy flux is in linear proportion to the trigger power. 

The natural unit ^ = c = 1 is used throughout in the present paper. 

II. DERIVATION OF QUANTUM AND SEMI-CLASSICAL EQUATION 

Consider three level atom (or molecule) of energies, e^ > eg > eg, as shown in Figdl]). We 
assume that transition between two lower levels, |e) and \g), are dipole forbidden. Suppose 
that the upper level \p) has substantial El rates both to \e) ,\g). (This can be replaced by 
weaker Ml transition, since the relation we need subsequently is the partial decay rate oc ef ■ 
with the energy level difference eij = e^ — ej, which holds both in El and Ml cases.) 



We focus on, and derive an effective hamiltonian of, two lower levels interacting with 
oscillating electric field E. Its hamiltonian density has been derived in |3|, [ij] for a single 
mode of field such as a light wave of definite frequency traveling in one direction. Extension 
to multi-mode fields such as counter-propagating modes of the same frequency is given in 
Appendix A. Its hamiltonian has a form of 2 x 2 matrix acting on two atomic states, |e) 
and \g), ~ E^AiE . The multi-mode field E may be decomposed into positive and negative 
frequency parts E^ = ^ . ^{E*e^^^^ + Eje~^'^^^) , where Ej, E* are slowly varying envelopes 
in time. We shall use variables, E^ = Eje~^^^* , E~ = E*e^'^^^ , to simplify formulas given 
below. In quantum field theory Ej and E* represent annihilation and creation operators of 
definite mode. The pertinent hamiltonian to our discussion of the single mode is 



(1) 



-ni= '^^ ''^' ' \ , (2) 

_ 2dpedpg _ 2dp^epa 

H'ge — ) H'aa — ^ _ 1 i V^ ~ \J i^j i \y) 

where |cep -|- jc^P = n{x) with n{x) the number density of atoms per a unit volume in a 
linear target region of < x < L. For simplicity we took isotropic medium and linearly 
polarized fields, taking E^ as scalar functions. The diagonal part oc ^aa of this hamiltonian 
describes AC Stark energy shifts, while off-diagonal parts oc /i^e are for two photon emission 
and absorption. 

For pll2 target the photon energy loq = eeg/2 ~ 0.26 eV is much smaller than level 
spacings to the electronically excited intermediate states, both Cpe and Cpg ~ 11 eV. Under 
this condition we may ignore uq compared to epa ,a = e,g in the formula for fiab and identify 




fiab to the polarizability for which precision calculation exists 18|. We thus use numerical 
values of parameters, figg ~ 0.80 ,/iee ~ 0.87,/ic,e ~ 0.055 all in the unit of 10~^^ cm^ [l9| 
for the pll2 Xv = 1 — )■ Xv = transition. 

The density matrix of pure atomic states, p = |^(a;, t))(^(x, t)| {{i/j{x,t)\ = (cg, Cj,)), obeys 
the evolution equation, dtp = —«[?//, p] . This quantum mechanical equation is generalized 
to include dissipation or relaxation. The needed variable, the density matrix for the mixed 



state, is given by a statistical mixture of pure states: 

p{x,t) = '^Ci\tpi{x,t)){4ji{x,t)\, ^Ci = l, 0<Ci<l, (4) 

i i 

with \ilji{x,t)) a set of ort honor mal pure state vectors. Dissipation occurs when a subsystem 
of |e) , \g) interacts with a reservoir and one integrates out reservoir variables due to our 
basic ignorance of the reservoir. The gen eral form of mixed state evolution including dissi- 



pation has been derived by Lindblad J20|], assuming the general principle of positivity and 
conservation of probability. As its result the time evolution equation of the density matrix 
has additional operator term, L[p\. The new additional dissipation term in the two level 
atomic system turns out equivalent to phenomenological relaxation terms given by two time 



constants, Ti,T2 (with the constraint Ti > T2/2 from consistency with 20]). 

It is convenient to write the evolution equations in terms of components of the Bloch 
vector defined hj R = ti pa = {iplaltp) . The basic Bloch equation including relaxation 
terms is 

dtRi = {pee - figg)E+E-R2 - tfige{e''^^'E+E+ - e-'''^^'E-E-)R, - ^ , (5) 

-'2 

dtR2 = -{Pee - Pgg)E+E-R, + /x,e(e"-*^+^+ + c-'^^^' E' E')R^ - ^ , (6) 

-'2 

dtRz = pge {i{e''''^E+E+-e~''^^^E-E-)Ri~{e'''^^E-^E~^ + e-'''^^^E-E')R2)-^^-^.{7) 

Ti ^ T2 usually, and the phase decoherence time T2 is much smaller and more important 
than the decay time Ti, which may be taken infinitely large for our practical purpose. 

Derivation of quantum field equation follows a similar line of reasonings. To perform 
the derivative operation d^ as in the Maxwell equation, one needs to calculate the double 
commutator; 

d^E^ = -[H,[H,E^]], H= f d^xCHf + iipHi) , (8) 

with the field energy density T-Lf = {E"^ + B'^)/2. For convenience we add less dominant 
oscillating terms of field modes to E"^ and use the locally well behaved field E{x,t) in "H/. 
The fundamental commutation relation in the radiation gauge QED [Ey{f,t) , Bzlf ,t)] = 



id^S^if — r') [2l| is used for derivation of quantum field equation. The result is 

{df - V^)E^ = V^VE^ , (9) 

-VE^ = (i^^l±tEln + ^^^1^R,\ i+ + p^,e-^^\R, - rR,)E- . (10) 



This equation 22| along with the Bloch equations ([S])~ dZD is the basis of the following 



derivation of our master equation. 

Slowly Varying Envelope Approximation (SVEA) Fast oscillating terms do not 
contribute to global features of time and spatial evolution when one makes averaging over 
a few X time and spatial oscillation periods. We thus extract terms that persist over time 
periods of typical light oscillation of order l/w both in time and space. Envelope functions 
denoted by Eji , E^ ought to be amphtudes of right and left moving components of rapidly 
oscillating parts oc e~*'^*^*^^'\ 

The result of SVEA may be summarized using dimensionless units of spacetime coordi- 
nates C, , T and dimensionless fields eL,R given by 

\E P B 

{^,t) = {amX,amt) , am{uj) = -^nfige{uj) , |eL,i?J^ = — bB—^ ri = —!-. (11) 

The quantity l/ctm = '^/ifige^egn) gives a fundamental unit of target length and time scale of 
evolution. Since a functional relation am{oj) = ami^eg—^) holds, the propagation problem of 
trigger irradiation of pair frequencies, u and teg — oj, is described by the same dimensionless 
quantities of a common a^- Its value at a; = eec,/2 is ~ 14cm and ~ 0.5 ns for para-H2 of 
density lO^'^cm"'^. 

The most general fundamental equations including both non-trivial propagation and PSR 
effects are derived in Appendix A and given by the formulas flA35p ~ flA42p . It is useful 



to recall the physical meaning of coupling constants fiat in the interaction hamiltonian, in 
order to fully appreciate the following approximation in our numerical simulations. Consider 
the extended hamiltonian including both of counter-propagating modes given by eq. (IA14p 
in Appendix. We first note that annihilation (oj) and creation (a|) operators of photon 
modes are related to complex fields by E^ ~ ai^/uj jIV , E~ ~ a\^uj/2V where V is 
the quantization volume. The important equations are obtained after SVEA and given in 



Appendix A. They are written in terms of envelope functions: 

+/i,e {{Ri - iR2T^El + (i?i - iR^f+^E],)\ , (12) 



[d, - d^)E, = ^ ((f^^l±tiln + ^^^^^Rf^)E, + ^^^^j^Ri-^En 



+ix,e {{Ri - tR2pE'i, + (i?i - iR2f-^El)\ . 



(13) 



The right hand sides of these equations give effects, all in bulk medium, of forward scat- 
tering oc ^^2£+A^^ + M.e-M33 ^(o)^ backward scattering oc ^'"'~^^'" R^f\ RL- pair annihilation 
oc fJ.ge{Ri — iR2Y'^\ and RR-, LL-pair annihilation oc iJige{Ri — iR2)^^^ ■ (The pair creation 
amplitudes appear in conjugate equations to those above.) Quantities R\ ^e^^*'^^ as defined 
by eq. (IA3ip are what are called spatial grating in the literature. The backward scattering 
terms, and RR,- LL-pair annihilation and creation terms are important only in the presence 
of spatial grating of polarization. Neglect of spatial grating is thus equivalent to retaining 
forward scattering and RL-pair processes, and ignoring all other terms. In the simple bound- 
ary condition set up below in this work, the backward Bragg scattering is expected to be a 



minor effect, and is also neglected in most works of SR. We refer to 23|] on the backward 



Bragg scattering effect in usual SR, and for instance to [2J] on the backward scattering 
effect on SR in low-Q cavity experiments. In a more comprehensive simulation in future we 
wish to quantitatively compute effects of the backward scattering, RR-, LL-pair processes, 
because non-negligible differences of these effects might arise in PSR unlike the SR case. 

In the rest of the present work we shall focus on PSR effects and ignore propagation 
effects which are much discussed in [Sj, [lj| and summarized in Appendix A. Explosive PSR 
events discussed below are expected to be insensitive to neglected propagation effects. The 
resulting Maxwell-Bloch equation for the single mode is 

drn = 47_(|efip + \eL?)r2 + ^'^{eReL)r^ - - , (14) 

drr2 = -47-(|eR|' + leiHn + 83fJ(e^ei)r3 - - , (15) 

a.rs = -8 (3fJ(e^ei)r2 + 9(e^ei)ri) - ^^^ , (16) 



1 2 

{dr + d^)eR = -(7+ + 7_r3)eij + -(ri - ir2)el , (17) 

(9^ - 9g)eL = 2 (7+ + 7-^3)eL + 2 ('^i ~ ^^2)e;^ , (18) 

7± = ^ • (19) 

Here Tj = amTi are relaxation times in the dimensionless unit. 

The dimensionless master equation fll4p ~ flTSl) is governed by two important parameters, 
the most important is T2 = amT2 and the next important is 7±. Another experimentally 
important parameter is the overall length and time 1/am oc 1/n, inversely scaling with the 
number density n. For larger number densities of excited atoms a smaller size target and a 
shorter time measurement of 0[ns] becomes possible. 

In terms of two component field if defined below the equation reads as 

{dr + a3di:)ip = -{-f+ + -f^r3)ip + -{ri - ir2)aiip* , y? = . (20) 

Magnitudes of R- and L-fiuxes change via RL mixing term such as 

{dr ± di:)\eR,L\^ = ri3(eReL) + r2$R(e«eL) . (21) 

R- or L-moving pulse alone propagates freely, because we ignored in this approximation 
non-trivial propagation effects. 

Quantum state of fields As usual in quantum field theory, we may interpret Er^l as 

annihilation and Ej^j^ as creation operator. The fact that the basic equation, (TT7|) ~ (ITS]) 

or (l20l) . simultaneously contains both annihilation and creation operators of field implies 

that the quantum state satisfying the field equation is given by a Bogoliubov transformation 

from the usual vacuum of zero photon state |0), 

00 
\^) = J2cn{x,t){EMr\0), (22) 

n=0 

where c„(x,t) is to be determined by (eq. fl^Ul) ) |\E') = 0. The quantum state |\E') is a mixture 
of infinitely many states of different photon number. We shall not pursue this line of thoughts 
any further, because we exploit the semi-classical approximation under the large quantum 
number limit of photons (the classical limit). The semi-classical equation is given by the 



10 



expectation value of quantum equation, 

{^\ Udr + (r,d^)v - ^(7+ + 7-r3)<^ " ^(^i " ^r2)a,^*^ |^) = , (23) 

with |\I/) the Bogohubov state given by eq. fl22p . The semi-classical equation turns out 
equivalent to replacing q-field operators in the quantum equation by corresponding c-number 
functions. Equations, flT^ ~ flTSj) . regarded as equations for c-number functions, thus 
constitute the master equation for polarization of medium and field. 

In our case of field condensate, medium polarization and fields are cooperatively involved: 
the target medium undergoes coherence oscillation, simultaneous with field oscillation, while 
keeping field envelopes slowly varying and finally almost time independent in a large time 
limit, as shown below. The field condensate part is technically equivalent to field state made 
of an infinite sum of multiple photon pair states in the so-called coherent state representation. 

III. IMPORTANCE OF INITIAL COHERENCE 

It is important to clarify in detail the ideal case of numerical solutions where all quantities 
in eq.( 1T^ ~ (TT8|) are of order unity, in the range of O[10 ~ 1/10]. For a deeper understanding 
of numerical outputs and a practical check of accuracy of numerical results, it is useful to 
know conservation laws of our non-linear system. We list in Appendix B all exact and 
approximate conservation laws that the system possesses. 

We have performed numerical simulations assuming CW (non-pulsed continuous wave) 
trigger laser irradiation of the same power from two target ends (called the symmetric trig- 
ger). This boundary condition is similar, but not identical, to the one of cavity mirror. Use 
of cavity mirrors has both advantage and complication. Two mirrors in cavity automatically 
generate counter-propagating waves, and they effectively increase the trigger power (which 
however is not critically needed in our case). On the other hand, each atom in cavity is 
affected by the same traveling fields many times and this complicates analysis. We use in 
the present work the simpler scheme of two CW counter-propagating triggers independently 
irradiated. 

Numerical results show the symmetric output fluxes from two ends, and we exhibit in 
the following figures one of these identical fluxes from one end. Result for the zero initial 
coherence ri(^, 0) = r2{C,, 0) = is shown in Fig([2]). Clear signature of delay much after T2 

11 



(~ 7T2 in this case) and explosive PSR is observed for strong trigger fields. It is difficult 
to obtain commercially available CW laser of this power. A reason of this difficulty is 
that relaxation of order T2 ~ 10ns may take over the coherence development under weak 
trigger usually exploited. Explosive PSR is a highly non-linear process having a definite 
trigger power threshold and disappears in this example certainly at the trigger power of 0.9 
MWmm~^, as shown in the inset of Fig(|2]). 
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FIG. 2: Time evolving output flux at a target end of length 30 cm resulting from the symmetric 
CW trigger irradiation of power 1 MWmm~ and 0.9 MWmm^ in the inset. (Note a large 
difference of ~ 10^ of the output power in two plots.) Assumed parameters are n = 1 x lO^^cm"^ 
of pH2, numerical values (see the text) of /i^fe for the transition Xv = 1—7- Xv = 0, relaxation 
times T2 = 10, Ti = 10^ ns's, and initial complete inversion (hence no coherence) of rs = 1 with 
''I) ''2 = taken for the initial target state. 



In Fig([2]) the complete inversion to the level |e) has been assumed as an initial condition, 
and it would be interesting to relax this condition and to further clarify neglected effects 
of the presence of initial coherence between two atomic levels, |e) and \g). There is an 
experimental method to imprint an initial coherence between |e) and \g) by adopting a 
clever excitation scheme. The atomic state right after excitation can be made a coherent 
mixture of two pure states, |e) and \g), namely Cp\e) + Cg\g) with |cep + |cgp = 1 at a single 



atomic site, by using the technique of STIRAP 25j. This kind of pure state may be formed 
by time overlapping excitation pulses of two frequencies, ~ e^e and ~ e^g. The state is a 
dark state, called so because no emission from \p) is observed despite of irradiation capable 



12 



of making both transitions, \p) — )■ |e) and \p) — )■ \g). 

The medium polarization r^ in the dark state is given by 



ri = 2a/p(1 -p)cos9o, ra = 2a/p(1 -p) sin 610, rg = 2p - 1 , (24) 

with p the fraction in the state |e). When this type of initial polarization of the dark state 
is formed, one may expect to expedite the coherence development for PSR, as shown in the 
following section. When CW laser is used for trigger, two overlapping pulses may induce 
PSR at the same time when the emission from \p) disappears: thus it may be called PSR 
from the dark. 

IV. NUMERICAL SOLUTIONS FOR HIGH DENSITY TARGET WITH INITIAL 
COHERENCE 

We first comment on what the number density n of target precisely means. This is the 
total number of atoms/molecules per a unit volume participating in PSR phenomena, hence 
it is the added sum of densities in the states, |e) and \g). Note also that the state \g) may 
or may not be the ground state of atoms or molecules. For instance, in the pHa transition 
of X^EJ"f = 2 — )■ 1, the target number density n may be much less than the ground state 
number density since \g) = {Xv = 1) is also an excited state. 

Time evolution from a dark state of initial polarization value given by eq. fl24p is illustrated 
for the pHa number density 1 x lO^^cm"^ in FigQ ~ Fig llTTl) . We exhibit dependence of 
the symmetric output pulse on the trigger power in the range of 10^^^ ~ lWmm~^ for 
n = 1 X lO^^cm^^ in Fig(j3j), which demonstrates two important features of explosive PSR 
with the presence of a large initial coherence: (1) the highest peak of PSR output is almost 
independent of the trigger power, suggesting a sudden, macroscopic release of energy (its 
density ^ ^eg^) stored between two levels, |e) and \g), (2) the onset time of explosive 
events, which may be called the delay time, depends on the input trigger power very weakly, 
and a linear logarithmic dependence has been confirmed up to 1 pW mm~^ (instantaneous 
enhancement factor ~ 8 x 10^^ in this case). A similar logarithmic power dependence of the 
delay time has been observed in numerical simulations of the single photon superradiance 
when the system is subjected to the trigger. 

The integrated flux is ~ |ii^maxpAt with At the time width of explosive event. l-EmaxP = 
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0[tegn] and this integrated flux is estimated as 0[l//ige], a quantity independent of tlie 
target number density ra, if tfie explosive event occurs. These figures show dramatic effects 
of initial coherence of the dark state. Observation of explosive events requires a target length 
> l/ttm oc 1/n. 

Detailed time structure of pulses as observed in Fig([3]) may differ if one adopts different 
available experimental parameters, but the output release of energy flux of order ee^n is 
universal in explosive PSR events. 

Spatial profiles of field fluxes and polarization components, Tj, within the target are 
illustrated in Fig(|l]) ~ Fig llTTI) . In this parameter set, about ~ 30 % of the stored energy 
eegU (the corresponding flux unit being 1.2 x 10^ Wmm~^(ra/10^^cm~^)) still remains in the 
target much after explosive PSR, and we observe a seemingly stable target state. Note that 
dimensionless fields Icjp = |-Ejp/(eegn) are plotted in Fig(jl]) and FigdS]). 
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FIG. 3: (Color online) Trigger power dependence of time-evolving output flux from the symmetric 
trigger irradiation of the power range, 10~^^ ~ IWmm^^, under the conditions of n = 1 x lO^^cm^'^, 
target length = 30cm, relaxation times T2 = 10, Ti = 10^ ns's, and the initial polarization, ri = 
l,r2 = r3 = 0. Depicted outputs from 1 Wmm^^ trigger power in solid black, from 10~^Wmm~^ 
in dashed red, and from lO^^^Wmm^^ in dotted blue are displaced almost equi-distantly in the 
first peak positions. Transition Xv = 1—7- Xv = of pH2 is considered. ~ 70 % stored energy in 
the initial metastable state is released in these cases. 

Result for the solid density of n = 2.6 x lO^^cm"^ and smaller excitation of r^ ^ —1 is 
shown in Fig( IT^ . There is a threshold of the excitation fraction of |e), located between 0.2 
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FIG. 4: (Color) Spacetime profile of dimen- 
sionless field energy, le^jp + le^p for the 1 
//Wmm~^ case of Fig ([3]). 




FIG. 5: Spatial profile at the latest time, 10 
ns after trigger irradiation, of FigQ. Note a 
large reduction by O[10^^'^] in the power scale 
in this figure. 
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1 /uWmm case of lig(13j). ^^ ' ^^— '^ 

% ~ 0.5 %, above which dramatic explosive PSR's emerge, as inferred from comparison of 
two plots of Fig([n]). 

So far we mostly showed explosive outputs in which most of the stored energy between 
|e) and \g) is released in a short time < 10ns after time delay. There is however a linear 
regime under a large initial coherence r « , -i = 1 , 2 in which the output flux is amplified in 
proportion to the trigger power. For instance, the amplification factor is ~ 10^ in the trigger 
power range of 1 fiWmm.^^ ~ 1 Wmm~^ for three different choices of initial r^ values of Fig 
(IT^ . In this figure we show the output fluxes in the linear regime taking as an example the 
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FIG. 8: (Color) Spacetime profile of r2 for the 
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FIG. 9: Spatial profile of r2 at the latest time, 
10 ns after trigger irradiation, of FigdS]). 
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FIG. 10: (Color) Spacetime profile of rg for ^IG. 11: Spatial profile of rg at the latest time, 

the 1 /iWmm-2 case of FigdS]). ^^ "^^ ^f*®'^ trigger irradiation, of Fig([l0]). 

trigger power of 1 mWinin"^. Although not shown in this figure, the hnearity of the output 
power to the trigger power has been checked for this set of parameters. 

V. STATIC REMNANT AND SPINORIAL SOLITONS 

In addition to dramatic explosive PSR emission it is also important to watch remnants 
after PSR emission, since previous figures at latest times may be taken to suggest formation 
of objects of non-trivial spatial profiles. Let us derive for this purpose the asymptotic form 
of fundamental equations. We anticipate that both the medium polarization f and fields 
cr, cl little change with time in the time region of t ^ 1/ttm after PSR emission. By taking 
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FIG. 12: Output flux for the solid target number density 2.6 x 10^^ cm^^ of length 2 cm, the trigger 
power 1 ;uWmm~^, relaxation times T2 = 10, Ti = 10"^ ns's, and smaller population r3 = —0.99 
(0.5% excitation), and r^, = —0.996 (0.2% excitation) in the inset. The other initial components 



are taken as ri = ^/\ 
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r2 = 0. Note a large flux scale difference ~ 10 in two plots. 
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FIG. 13: (Color online) Output flux of weak pH2 PSR in the linear regime in which the output 
power ~ O[10^] x the trigger power, for initial {r^, ^i, ^2) = (0, 1, 0) in solid black, (l/v2, l/v2, 0) 
in dashed red and (— 1/\/2, l/\/2)0) in dotted blue, using the same set of other parameters: n = 
1 X 10^'^cm^'^, target length = 1.5 m, relaxation times T2 = 10, Ti = 10^ ns's, and the trigger power 
1 mWmm~^. The output power scales with the trigger power, as explicitly checked in the range 
of 1 yuWmm"^ ~ 1 Wmm~^. 
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vanishing time derivatives, one may eliminate polarizations Tj in favor of field components 
and write profile equations of spatial variation for fields, 

e'ij = "iigcR + if el , e'^ = -2igeL - ife*^ , (25) 

g = g{eR, cl) = 7+ - 7_ _,„ ^ 2/1 — fTTI — \2\'2 , ra \ i2~rT ' (^b) 

167lr|(|eij|2 + \eL\^y + 64rir2|eijeLP + 1 

f ^ f, .^ ^riCRCL ((47_r2(|eHp + \eL?) - 2) , . 

^ ^^"^' "^^ 167V|(|e^|2 + |ez.P)2 + 64rir2|e^ez.P + 1 ' ^ '^ 

"■' 167^r|(|e^J2 + |eLP)2 + 64rir2|efieiP + l' ^ '' 

where ' indicates the spatial derivative d^. 

Despite of complicated field dependent coefficient functions that appear in /, g, the struc- 
ture of profile equation fl2^ is rather simple. Oscillatory behavior governed by terms oc g 
can be eliminated by taking three bilinear forms of fields, le/jp , le^p , eRCL. 

{\eR? + \eL?)' = , (|e^p - \eL?)' = -453(/e^eI) , (e^e^)' = -^/(|e^p - \eL?) , (29) 

where the function / depends effectively on crCl alone since the total flux is a constant of 
integration due to the first equation of fl2^ . hence with a real constant Cq, |ei?(0P+|6L(0P = 
Cq. The set of profile equations, fl29|) . is transformed into two equations of phase functions, 
V5(0>'S'(0, defined by 

ei?(0 = eo cos ^(0 , 6^(0 = eoe^^(«) sin <^(0 , (30) 

^' = 1^1, 2 4 2!l!i% r^T^ sin(2^) , (31) 

1 + 167iegr2^ + 1665x1X2 sm^(2v9) 

S' = ,^,,2 4f2fif -^7^ ^°<2^) ' (32) 

1 + 16716^X2^ + 166(5X1X2 sm^(2v9) 

with (p{l/2) = IT /A , S{l/2) = 0,1 = amL. Since 6j's contain four real functions, the result- 
ing two equations here reflects a non-trivial self-consistency of the ansatz fl30|) . A similar 
equation with R -H- L interchanged may be set up, suggesting another kind of solitons. 

Equation for the angle function ip{C,) fl?I]) is self-contained, and has the following analytic 
solution under the boundary condition 6ij(//2) = 60/v^: 



2 , 1 67^x|6g + ^ In "^ ^ ~ ^/^ 

32x1X26^ el - e\ 4xi 



24 - 4 + :;A"4 i^7r^ = -^^^r^ • (33) 
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Field may decrease exponentially in the central region, like e^ oc exp[— 8x260!^ — 
//2|/(167^r|eQ + 1)]. One may define the soliton size by the e-folding factor as ^^ = 
(167^r|eo + l)/(8r2eo). The actual soliton size is Xs = ^s/c(m- 

The spatial variation of cr oc cos(f{C,) is monotonic, decreasing or increasing depending 
on the if region of either [0,7r/2] or [tt/2,tt] (defined modulo vr). These two fundamental 
regions are separated since ip' = at edges of these regions due to sm{2ip) = there. One 
may identify these two solutions as different objects. Either of fields cr , e^ vanishes at edges 
of fundamental regions, but not both. Solution defined by the fundamental region [0, 7r/2] 
corresponds to absorber soliton in which both R- and L-fiuxes are absorbed at edges, but not 
emitted at the other edges, as illustrated in Fig flT^ and Fig( lT5|l . This object may be called 
absorber soliton. The other fundamental region [vr/2, vr] gives emitter soliton which may be 
realized when the excited |e) state is sufficiently occupied. The existence of two types of 
soliton condensates is an important result indicating existence of a new kind of topological 
soliton whose topology is discussed in Appendix C. 

When the target size L is large and L ^ C^s/ctm, one may expect a copious production 
of absorbers and emitters within the target. When the target size is smaller than ^s/otm, 
the target edge effect becomes important (in general destroying, or blocking its formation 
of, soliton), and it may be difficult to create a soliton. 

400000 r 
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FIG. 14: (Color online) Profile of fields and r^ of helical absorber soliton. |i?_Rp (in red), \Ei\^ (in 
dotted blue), both in the unit Wnini^^, and r^ in arbitrary unit (in dashed black) are plotted for 
a case of n = 2.6 x 10^^cni~^,T2 = 20ns, Ti = lO'^ns. r^ ~ —1 near edges and ra w —0.8 in the 
middle. 
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Soliton solution obtained by direct numerical integration of (IHT]) is illustrated in Fig( lT^ 
along with distribution of the population difference r^. Solitons are characterized by two 
end points of rs ^ —1 and an intermediate region oi r^ ^ 0. It is important to have a long 
enough target in order to accommodate many solitons within the target. Soliton size can be 
made smaller if one can use a larger target number density close to the solid density. 

VI. CONCLUSION 

In summary, we derived and numerically solved the master equation for time evolution 
of PSR emission and formation of field condensates in long dense targets. We have demon- 
strated (1) numerical identification of two different types of PSR events, explosive and weak 
ones, and (2) theoretical existence of spinorial solitons stable against PSR emission. Realistic 
experiments can be designed using numerical solutions of our master equation. 

Note added. 



Recently, we became aware of a related work 26| where the time evolution of triggered 
two-photon coherence is examined. The authors of (26j treat the field differently from the 
one of our semi-classical approach, which results in our coherence development time of order 
several nano seconds in dense targets, much shorter than their value. 
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Appendix A: Two level effective model interacting with multi-mode fields 

We extend results of Appendix A in [ij] to the case of multi-mode fields such that two 
color problem including propagation effect is properly treated. This is the most general case 
of two photon problem. Its notation in this reference is slightly changed. 

Atomic system 
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The state vector of an atom can be expanded in terms of the wave function, 

|^(t)) = c,(t)e-«*|(7) + Ce(t)e-=*|e) + ^c,(t)e-^» . (Al) 

p 

Ca(t)^s are probabihty amphtudes in an interaction picture where e^'s are energies of atomic 
states. 

The atomic system may interact with hght fields. The electric field E{x,t) that appears 
in the hamiltonian via El or Ml transition is assumed to have one vector component alone, 
namely we ignore effects of field polarization. This is a valid approach under a number of 
circumstances. One then decomposes the real field variable E{x,t) into Fourier series, e"*"^^* 
times a complex envelope amplitude Ej{x,t), and its conjugate, where Ej{x,t) is assumed 
slowly varying in time, 

E{x, t) = J2 (^j (^' ty'^'' + Eji^^ t)e-^^^*) . (A2) 

i 

Each discrete mode j is taken independent. The most interesting are the cases of two modes 
with Ui + UJ2 = eeg and the single mode with u = eeg/2. 
The Schrodinger equation for a single atom, 

t^^\m) = iHo + dE)m)), (A3) 

with Hq the atomic hamiltonian, is used to derive the upper level amplitude Cp{t). Using 

t^^{p\m) = {p\{Ho + dE)mt)), (A4) 



one has 



.dcp 



i^e-''"' = {dpeCee-''^' + dpgCge-''^')E , (A5) 



where dab are dipole matrix elements. This can formally be integrated to 
Cp{t) = -i dt' (dp^Ce{t')e''^^*' + dpgCg{t')e''^^^'\ E{x,t') 

dt'J2 {dpeCe{t')e''^^'' + dpgCg{t')e''^'''\ (E*{x,t')e'''^'' + E,{x,t')e^"^^'') ,(A6) 



with eab = €a — eb the atomic level energy difference. The initial condition Cp(0) = is 
assumed here. 
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Markovian approximation and effective two level model 

The basic strategy of deriving equations for tlie lower two level amplitudes Ce, Cg in a closed 
form is to eliminate atomic variables Cp's related to the upper levels. This is essentially done 
by neglecting a long-time memory effect (the Markovian approximation) and making slowly 
varying envelope approximation (SVEA). The idea of the Markovian approximation is to 
replace dynamical variables, Ce(t'), Cg{t'), Ej{x, t') in the integrand of eq. (IA6p . by their values 
at time t, neglecting all the past memory effects. This gives 

Ef \ — Qi('^j+<^pe)t ^ _ ^-i(u}j-epe)t 
dpeCe ■ E* Ej 



+dpgcg { —-— E* - ---— E, ] , (AT) 



which is inserted into equations for the lower levels 



-^ = -tJ2 depEix, t)c,(t)e-^^-* , (A8) 



dt 

dCg 



dr 

^ = -^ E ^5P^(^' ^)cp We-^^-* . (A9) 

p 

Note that dab = dba are real by an appropriate choice of phases. 

We ignore rapidly oscillating terms, keeping in mind the two most important cases of the 
mode choice. The result is 

d Ce 



Cg 




,,. , ^-^.. ., (AlO) 






n 



SjMjj'Sf , 



. ^ ^^ dl^{2epe + ujj-ujf) ^^ dl{2epg + Uj - ujj>) 

fieeiUj, U,>) = 2^ — - , figgiu^U,,) = }_^ — , (A12) 

p ICpe + ^j)[(^pe - ^j') p [^pg + ^j)[^pg " ^j' ) 

s^ dpjpg{epg - {uj + ujj')/2) s^ dpjpg{epe + {ujj + ujj')/2) 

lieg{UJj, UJj') - 2_^ 7- -TT; — ^ , f^ge[UJj, ^j<) - 2_^ , i , , V^ i ,, ^ • 

^ iCpg - UJj )[epg-ujf) ^ [epe + UJj )[epe + ujj> ) 

(A13) 

S may contain both Ej and E*. We simplify notations below such that fields are redefined 

incorporating oscillating factors e"*'^^* in E^ etc. 
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Single color problem 

We apply the result to the problem of counter-propagating fields Eji and Ei of a single 
color of Uo = eeg/2. In the 2x2 hamiltonian (lAllI) one may use the complex field E^ = 
Eje~^'^^^ (the positive energy component corresponding to the photon annihilation operator) 
and its conjugate E~ = E*e^'^^^ (the negative energy component corresponding to the photon 
creation operator) to eliminate phase factors except e^*^'=»*, as is done in the main text. Since 
each mode is independent, it separately satisfies the field commutation relation, necessary 
for derivation of the quantum field equation, justifying the result of manuscript for the 
degenerate case. 

More concretely, 

_^ , \ee{ElE], + EtEl + EtEl + EtE],) e--*/i,e(i?A^A + ^l ^l + '^E^rEl 




e— V.e(i?^i?^ + ElEl + 2E'^E-^) f^gg{E^E], + ElEl + E^El + ElE'^^ 

(AM) 

_ ^d.ped'pg _ ^"'pa^pa (' A 1 c;^ 

^pg + £pe ^pa " ^0 

(RR) and (LL) terms describe pulse propagation with compression and splitting, while (RL) 
terms back-scattering, pair creation, and pair annihilation. 

Two color problem 

We may consider for Er , El envelopes of two different colors of uji with ui + uj2 = Ceg- 
Separation of cross mode terms leads to 

'Hi = T-Ld,R + 'Hd,L + 'Hi2,R + 'Hi2,l + 'Hi2,rl , (A16) 

^_| Ei,af^ee{0Ja,UJa)E+,E-^ e'^'"' Ei,a f^eg{oJa, (^a)E+,E+ 

e~'''''Ei,af^9e{^a,UJa)E-^E-^ Ei,a /^^^(^a, Wa)^,>^,> 

_' fXee{uJi,U2)E^^E-^+ fXeeiuJ2,(^l)E^2E-^ e^^^«*2/ieg(Wl, ^2)^+1^,^2 

-rii2,i — I _ _ _ 

e-*^^«*2/i3e(wi, a;2)^i,i^i,2 f^ggi^i^ ^2)E^^E-^2 + /i<,<,(w2, Wi)S+2^i,i^ 

(A18) 

i'Hl2,RL)aa = ^^aa{^l■,^2)E ^^^Ej^^ + ^^aa{^2■,^l)E ^2E L,l 

+/^aa(Wl, U2)EJ;^^E]^2 + flaa{(^2, UJi)E^.,E]^^^ , (Al9) 

- {'Hl2,RL)eg = e^^-*2/Xe5(^l, C^s) (^^,i^^,2 + i?A.2<i) , (A20) 
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dla{2epa + uji- UJ2) 

/iaa(Wl, OJ2) = 7 ■ rj T , (A21) 



l^gei^, Ceg - w) = Hgei^eg - U , u) = fieg{u}, e^g - Uj) = fieg{(^eg " W, Uj) 



Bloch equation 

The Bloch vector defined by 



^pe^pg\^pg ~r ^pe) 
[tpe + Uj){tpg -U) ' 

(A22) 



'^e'^e C^Cg 



i?=(^|a|^)=tr<Tp, p=|V;)(V^|= ^^ ^M, (A23) 



'"s'"^ '^g'^9 



satisfies quantum mechanical equation (disregarding relaxation terms) dtR = —itr cr[Hi , p] . 
Explicit calculation using the Hamiltonian above gives 

dtRi = (/iee - figg)E+E-R2 - ipge{e''^^'E+E+ - e-''^^'E-E-)Rs , (A24) 

dtR2 = -(/iee - Pgg)E+E-R^ + /i,e(e*^-*^+E+ + 6-^'^^' E' E-)R^ , (A25) 

dtR^ = fige {i{e''^'^E+E^ - e~''^^^E~E~)Ri - (e^^=«*E+E+ + e'^^^'^E- E-)R2) .(A26) 

We suppressed mode index j for simplicity. Conservation law holds: dt{R\ + R\ + R\) = . 

Field equation 

Commutation relation of fields necessary for derivation of quantum field equation 
[Ey{f,t) ,Bz{f',t)] = idxS^{f— f') , is valid for each independent mode. The double com- 
mutator, 

d^E^ = -[H,[H,E^]], H = I SxiUf + tr pUi) , (A27) 

il pHl = {^l^\H\t^) = -(/ieelCeP + Pgg\Cg\^)E+ E' - Pge{clCgE+ E+ + c^CeE- E') , (A28) 

with the field energy density "H/ = {E"^ + -B^)/2, is calculated as 

{dl-V^)Ef = V\V,,,E,,)^, (A29) 

-^,T^; = (^J^^^±^M^n + il^^l^MlLR^ El, + (/i,,),,,e--*(i?i -^i?2)^^ • (A30) 

SVEA and dimensionless equations for two color modes 

All terms both in the Bloch and field equations must have the same oscillatory behavior 
for global evolution of polarization and fields. This gives a phase matching condition of the 

24 



form Ui + U2 = teg and momentum balance with E^ oc e**^^ , E^ oc e"*"^^. For time SVEA 
one may then ehminate the phase factor e^*^'^^* in the Bloch equation. For space SVEA we 
introduce spatial variation of polarization of the form, 

R^ = i?f + i^j+^e'^"" + i?;"^e-2^'^" . (A31) 

LHS of field equations ~ —2iuj{dt ± dx)Eji^L for the counter-propagating modes of the 
same frequency, hence (with d± = dt ± d^) 

d^En = Y (( ^"^^"" n + ^''~^'' Rf^)En + ^^^^^J^R^^El 

+/i,e {{Ri - iR2T^El + (i?i - iR2i-'^E*^)\ , (A32) 

d_EL = — f^i^SlU^n + ^''~^'' Rf)EL + ^'^'^'' R^i^En 

^ \ ^ Zi Zi 

+/i,e ((i?l - iR2TEI + (i?i - 2i?2)(-)^l)') . (A33) 



We introduce the dimensionless unit: 



e 



^(l),(2)|2 p 



(^,r) = (a„a;,a„t), a„,{uj) = '^n^ge{uJ,teg-u) , Wl^r ? = — — — '- , r^ = — . (A34) 

Assume R-mover of frequency Ui and L-mover of frequency U2 (neither R-mover of frequency 
UJ2 nor L-mover of frequency Ui). Note the universal parameter HegioJi, W2) = l^geij^i, ^2) for 
any combination of uji + uj2 = teg- The master equations for medium polarization and fields 
are 

(0) 

-2z(ei - (6^)^)4+) - nel - {elf)r{'^ - ^ , (A35) 



T2 



9.rW = 47i^^^e«elrf -2^(e?,-(el)2)rf +4(7i^)|e«P + 7i'^|eLr)4+^ + 8$>(eKe,)rW "^ 



(+) 



r2 
(A36) 

A (0) /l^ (1)1 |2 I (2)1 \2\ (0) , oWi^ \ (0) A (12) * (-) A (21) * (+) 

^(0) 

+2(ei + (e^)^)rW + 2(e| + (el)^)^-) - ^ , (A37) 



^2 

5.rW = -47i^^)e^elrf)+2(e?,+ (el)2)rf -4(7«|e^P+7i'VLnri+^+85R(e^e^)r(+)-^ , 

T"2 

(A38) 
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9,rf = -8 (3fJ(e«ei)rf + 53(e^e^)rS°)) + 2z{el - (eDyt^ + 2%{el 



2V(+) 



^R) }'l 



SO) 
-2{ei + {e^^fyf' - 2(e^ + {eD'Yp - '-^-^ , (A39) 






a.r5+) = 2zrS°)(e^-(e2)2)-2rf(e^+(eI)2)-8(3fJ(e^e^)r(+) + 53(e^ei)ri+))-^, (A40) 



(+) 
3 



(5. - d,)e, = ^(7f + 7i^Vf )e, + '-^'r^'e^ + ^(r^ " -f )e^ + ^(^i"^ - -r))el . 



(A41) 

M- 
(A42) 

(a) _ /iee(l^a, CUg) ± l^ggj^a, ^a) {ah) _ /iee(l^a, l^b) ± l^gg{.^a^ ^h) ( ^A'X'\ 

^f^ge ^fJ'ge 

Oi = , aij = — ^ , (A44) 

^eg ^i^eg 

with yUdb defined by (lAlSp . 

The single mode equations in the text are readily derived by taking Ui = eeg/2, ai = 
1, aij = 1 and all 7^^ a, 6- independent. 

Pulse compression factor 

We shall estimate pulse propagation effects neglected in the text. Pulse propagation may 
be described by ignoring RL mixing terms in the general master equations. By taking one 
mode cr of one color, the basic propagation equations are 

9,ri = 4r33e^ + 47_r2|eij|2--, (A45) 

drr2 = Ar^?fiel - A^_n\en\^ - ^ , (A46) 

T2 

d^r, = -Ain^el + r^^el) - ^-^^ , (A47) 

{dr + d^)eR = '- ((7+ + 7-r-3)eij + {n - ir2)e*ji) . (A48) 

We shall ignore relaxation terms, taking r^ — > 00. Results of [iJ] in terms of the area 
function follow with the assumption of reality of the function cr. The relation ri = — 7_r3 
automatically follows from the consistency of three Bloch equations. The fundamental equa- 
tion of the propagation problem is given by a single non-linear field equation in terms of the 
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area function 9{^, r): 

9 drO cos 9 ^ / . X 

^«=77t^=^' ^3 = ± /— ^ , r2 = ±sin0, (A49) 

(9, + (95)(9,^ = ± sin edrO . (A50) 

Analytic solutions of this non-linear equation give [ij] 

(1) Pulse splitting. The number A^ of split pulses is given by the pulse area of the initial 
flux Fi{t) divided by 2n: 

N = ^^f^je + if^ee-f^jy^l^ dyF,{y) . (A51) 

(2) Pulse compression. The pulse of area < 27t is compressed by an amount E (result 
obtained for Lorentzian pulse), 

E = . . —. — , (A52) 

{amx sin(^/2) ± cos(e/2))2 + sin2(^/2) 

9 = J/x2^ + (/Xee-/x,,)V4 / dyF,{y) , (A53) 

J —oo 

± depending on amplifier (absorber). 

We may estimate the pulse compression factor flA52p for CW trigger irradiation of dura- 
tion t in which case 6 ~ /3t: 

E = ^ ~ ^ . (A54) 

{PtamX ± ly + {(3ty l±2l3ta^x ^ ' 

In all cases of our interest [it < [3T2 ^ 1. Thus, unless the target length is large enough, 
close to l/(2/3T2am), the effect of pulse compression is not large. 



Appendix B: Exact and approximate conservation laws 

We focus on the degenerate case of ui = U2 = ^eg/2. There are three different classes of 
exact and approximate conservation laws: (1) one exact conservation that holds with finite 
Ti, (2) one more approximate conservation law that holds in the Ti — > oo limit, (3) one 
further approximate conservation law that holds in the T2 — )■ oo limit (Ti ^ T2 assumed). 

The first exact conservation law is derived directly from two equations of motion for the 
field Cj and it reads as 

{d, + d^)\eR\^ = {dr-d^)\eL?. (Bl) 
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An integral form of this conservation for a finite target of length L {I = UmL below) is 

^ f^ d^M' - \eL?) = -[\eB\' + le^l^U • (^2) 

The integrated quantity of |eijp — |eLp stored in the target balances against its flux outgoing 
from two target ends. For the symmetric trigger, RHS of this equation vanishes, and the 
integral in LHS is a constant of motion. 

The second conservation law that holds in the Ti — ;■ oo limit is 

dr{n + 4(|efi|2 + le^p) ) + 49e(|e«|2 - le^l^) = , (B3) 

corresponding to the energy conservation. The energy density inside the target is a sum of 
medium and field energies, r3/2 + 2(|eijp + le^p), in our dimensionless unit. The integrated 
form of this conservation law in the real unit is 

d_ '■^ 
dt 

The third class of conservation law that holds in the T2 — ?■ oo limit is 



j dx {^Rs + 2{\En\' + \El\')) = -2[\En\' - \EL\']to ■ (B4) 



drirl + rl + rl)=0. (B5) 



Appendix C: Helical soliton 



A new type of topological solitons may exist, because the basic equation has two compo- 
nents v^^*^(0 j"^ = 1;2, and one can give a topological quantum number in 1+1 dimensions, 
as illustrated in Fig flT5|l . For simplicity assume two real component field {X,Y) and its 
periodicity with period of the t arg et length /(= amL) or a few times of this length. We 



may define the homotopy class 27|] of the mapping of a circle x + iy = /e*^ , < .^ < 27r in 
two dimensional real space onto the field space of the unit magnitude, X"^ + ¥"^ = 1. The 
winding number w is defined using the complex field Z = X + iY = e*"^*-^-*: 

>'2tt 



W = —I 



r ^z-a,z = ^<2^1^^ . (ci) 



When this winding number is quantized, w = n ,n = 0, ±1, ±2, ■ ■ ■ , the winding number is 
topologically stable and conserved during time evolution. 

The correspondence to static solutions in Section |V] is as follows. One considers the 
real 3- vector field X{C,) of unit length, {X,Y,Z) = (cos v?, cos 5* sin y?, sin 5 sin y?) with {p,S 
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FIG. 15: (Color online) Helical structure of absorber soliton. Target region is irradiated from 
both ends by trigger lasers of different colors. The X(^) = (cos (^, cos 5 sin 99, sin 5 sin 99) with ip^S 
defined by eq. ()30p may wind. In the return trip from the right edge to the left edge, not shown 
here, X{^) further winds and comes back with ^(vr) = — X(0) at the left edge, giving a spinor 
field. This is an absorber soliton without emission at two ends. 



identified as the phase variables in static solutions, and a mapping of unit circle < .^ < 27r 
onto X{^) space. Two solutions of eg. fl3T]) . f l32|) corresponding to two different fundamental 
regions, [0, 7r/2] and [7r/2,7r], are connected together at ^ = ti/2. Then, in the return trip 
of .^ = 7r/2 — )■ Tx from the right edge to the left edge of soliton, the orientation of X{S,) 
is further advanced forward (dictated by continuity of solutions), and finally comes back 
with X{'k) = — X(0) at the left edge. This means that solutions of X are two-valued 
representation, namely spinors. 
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